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We introduce a large class of holomorphic quantum states by choosing their normalization func- 
tions to be given by generalized hypergeometric functions. We call them generalized hypergeometric 
states in general, and generalized hypergeometric coherent states in particular, if they allow a resolu- 
tion of unity. Depending on the domain of convergence of the generalized hypergeometric functions, 
we distinguish generalized hypergeometric states on the plane, the open unit disk, and the unit 
circle. All states are eigenstates of suitably defined lowering operators. We then study their photon 
number statistics and phase properties as revealed by the Husimi and Pegg-Barnett phase distri- 
butions. On the basis of the generalized hypergeometric coherent states we introduce generalized 
hypergeometric Husimi distributions and corresponding phase distributions as well as new analytic 
representations of arbitrary quantum states in Bargmann and Hardy spaces. 
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1. INTRODUCTION 

Coherent states, their variants and generalizations 
have been extensively studied over the last four decades. 
A comprehensive review of this development can be found 
in Refs. Recently, Klaudcr et al. Q have ex- 

posed a rather general method for constructing holomor- 
phic coherent states of the form 



1^) = m\^\')r'^ E 



(1.1) 



where z is a complex variable, p{n) a set of strictly posi- 
tive parameters and the states |n) form an orthonormal 
basis (usually the Fock basis). The normalization func- 
tion is given by the scries 



— ' p(n) 



(1.2) 



for C = its radius of convergence determines the 

domain of definition of the states in Eq. (|l.l|l . 

In this paper we consider a particular, but neverthe- 
less very large class of holomorphic coherent states of the 
form Hl.l|l . for which the normalization functions (|1.2() 
are given by generalized hypergeometric functions 



■^(C) = pPqk^^-. ■■■,a,p;bi,-- • ,&g;C)- 



(1.3) 



We call them generalized hypergeometric states in gen- 
eral, and generalized hypergeometric coherent states if in 
addition they allow a resolution of unity with a positive 
weight function. They comprise many well known quan- 
tum states as particular or limiting cases as well as some 
of the examples constructed in Ref. 3 . 



The plan of the paper is as follows. In Section |21 we 
introduce the generalized hypergeometric states, which 
naturally fall into three classes: generalized hypergeo- 
metric states on the plane, the open unit disk, and the 
unit circle. In Scction|31we discuss the generalized hyper- 
geometric coherent states and use them to represent an 
arbitrary state by analytic functions cither on the plane 
or the unit disk. We then show in Section 0] that all 
states are eigenstates of suitably defined lowering op- 
erators. In Section \S[ we consider the photon number 
statistics of the generalized hypergeometric states and in 
Section |S1 their phase properties as given by the Husimi 
and Pegg-Barnett phase distributions. Generalizations 
of the Husimi distribution on the basis of the generalized 
hypergeometric coherent states and corresponding phase 
distributions are introduced and discussed in Section |7| 
We conclude in Section |H1 with a summary of our main 
results. 



2. GENERALIZED HYPERGEOMETRIC 
STATES 

In this section we introduce the generalized hypergeo- 
metric states, classify them according to definite criteria 
of convergence, and indicate the constraints on the pa- 
rameters on which the states depend. 

We define generalized hypergeometric states (GHS) by 



\p;q;z) 



\ai, 



, flp, &l, . . . , bq, z)(^p.q^ 



[pAA,(|z|2)]-^ £ 



\/pPq{'^) 



|n), (2.1) 



with the strictly positive parameter functions (of the dis- 
crete variable n) 
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pPq{n) = ppq{ai, . . . ,ap;bi, . . . ,bq;n) 

{bl)n ■ ■ ■ {bq)n 



T{n + 1) 



{ai)n ■ ■ ■ {ap)n 



(2.2) 
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where (a)„ = T{a + n)/T{a) is the Pochhammer symbol. 
The normahzation function is then given by the general- 
ized hypergeometric function 



P-^qiO = pFq{ai, . . . ,ap]bi, . . . ,bg;C) 
(ai)„ • • • (ap)„ C" 



n— 



(2.3) 



evaluated at ^ = The GHS depend on the complex 
variable z and on the two sets of numerator (ai, . . . , Op) 
and denominator {bi,...,bq) parameters. The non- 
negative integers p and q assume the value p = (g = 0) if 
there are no numerator (denominator) parameters; this 
will be indicated by a dot, e.g. |l;0;z) = |a; • ; z)(x.o). 
Since the parameter functions pPqin) arc the defining 
quantities, some of their particular properties become 
"global" and hold true for the GHS and all quantities de- 
rived on their basis. This will by exemplified repeatedly 
in the following. As a first property we mention the sym- 
metric dependence on the numerator and, separately, the 
denominator parameters. A second property is related 
with the factorized structure of pPq{n): if / numerator 
and / denominator parameters coalesce {at ~ bj), they 
cancel and the state \p; q; z) reduces to the lower-order 
state \p — I; q — I] z) of the remaining parameters. Thus, 
the conventional coherent states (CS) result for p = q = 
or, more generally, for p = q ^ and all parameters co- 
alescing. 

The normalization function also determines the scalar 
product of two GHS with identical parameter sets 



{p; q;z\p;q;z) 



pAfq{z*z') 



^pNq{\z\^) pNq{\z'\ 



(2.4) 



It follows that the GHS are normalized, but not orthogo- 
nal. The scalar product is well defined if the generalized 
hypergeometric functions involved converge. The func- 
tion pFq[ai, . . . , ftp; 6i, . . . , bq] C) converges in the follow- 
ing cases for: 



any finite C, 

ICI<i 
ICI-1, C^i 



if p<q + l] (2.5a) 

if p = q+l] (2.5b) 

if p = q + l, ri<Q; (2.5c) 

if p = g + 1, < 77 < 1, (2.5d) 



where 77 = Re(^ 



3=1 «J 



_ Sj=i^j); i'^ other cases 

it diverges According to the three cases of uncondi- 
tional convergence, Eqs. H2.5a|l to (|2.5c(l . we obtain gener- 
alized hypergeometric states on the plane, the open unit 
disk, and the unit circle, respectively. These states de- 
pend analytically on the complex variable z (disregarding 
the normalization function) and thus belong to the very 
large class of holomorphic quantum states. Besides these 
normalized states we shall consider also unnormalizable 
states on the unit circle violating the condition 77 < 0. 

We conclude this section by reviewing the constraints 
on the parameters ai, . . . , ap and 61, . . . , 5^. First of all, 
we exclude zero or negative integer values for ai, . . . , Op 



(leading to undefined parameter functions pPq) and for 
bi, . . . ,bq (leading to undefined normalization functions 
pFq). Next we require the parameter functions pPq{n) to 
be real and strictly positive, leading to the conditions 



jbi+n)--- [bq + n) 
{ai+n) - ■ ■ {ap + n) 



> for aU n = 0, 1,2, 



which follow from the obvious recurrence relation 

pPqjn + 1) _ pPqjn) (bi+n) - ■ ■ {bq + n) 
n\ (fli + n) ■ ■ ■ {ap + n) 



1)! 



(2.6) 



(2.7) 



with seed value pPq{0) = 1. The conditions (|2.6|) are sat- 
isfied for all n, if: (i) all parameters are real and positive, 
(ii) an even number of parameters is real and negative 
with pairwise the same negative integer parts, (iii) an 
even number of numerator parameters and/or of denom- 
inator parameters are pairwise complex conjugate to each 
other. If cases (ii) and/or (iii) occur, the remaining pa- 
rameters are supposed to be real and positive. 



3. GENERALIZED HYPERGEOMETRIC 
COHERENT STATES 

In this section we introduce generalized hypergeomet- 
ric coherent states (GHCS) as a set of generalized hy- 
pergeometric states (GHS) which is (over) complete and 
allows a resolution of unity with a non-negative weight 
function. The resolution of unity is then used on the 
one hand to determine the corresponding weight func- 
tions and on the other hand to introduce new analytic 
representations for arbitrary quantum states. 

We consider first the GHCS on the plane and on the 
unit disk related, respectively, with the convergence cri- 
teria (|2.5a|) and Ij2.5bp . In these cases the resolution of 
unity (1) takes the form 

i j d^z pWq{\z\'^) \p]q]z) {p]q]z\ = i, (3.1) 

with a positive weight function pWq{x) depending on 
X = the integration is over the complex plane or 
unit disk, as appropriate. Introducing the states H2.1|) 
into Eq. H3.1|) and performing the angular integration, 
the following conditions result: 



dx pWq{x) = pPq{n), 71 = 0,1,2, 



(3.2) 



where pWq{x) = pWq{x) / pMq{x) and the upper limit of 
the integral is i? = oo {R = 1) for the GHCS on the 
plane (disk) . Correspondingly, the unknown distribution 
pWq{x) is given by the solution of a Stieltjes (Hausdorff) 
moment problem with the moments given in Eq. (|2.2|l . 
As shown by Klauder et al. the solution can be ob- 
tained by using Mellin transform techniques. Thus, re- 
placing in Eq. H3.2|l the discrete variable n by the complex 



3 



variable (s — 1), the distribution pWq{x) and the parame- 
ter function pPq{s — 1) become a Melhn transform related 
pair. There are well known references where such pairs 
are tabulated [IIIH. Thus, using Ref. H p. 303, for- 
mula (37)] or Refnir P- 728, formula (8.4.51.9)], we find 
in terms of the Meijer G-function 



R 



dx 



-1, 



^p, q+l I ^ 



ai - 1, 



r(.) 



r(s + 6i - 1) 



. . . , ap - 1 
...,6,-1,0 

r(s + 6,- 1) 



T{s + ai-l)---V{s + ap-l) 



where the r.h.s. is the s-dependend part of pPg(s— 1). It 
follows that the weight function in Eq. H3.1|) is given by 



pWq{x) = pWq{ai, . . . ,ap;bi, . . . ,bq;x) 

_ r(ai) . . .r(ap) 



Tih)---T{bq) 



pFq{0'l, ■ ■ ■ , O-p] bi, 



fli - 1, . . . ,ap - 1 
6i-l,...,6, -1,0 



The Mellin transform (|3.3|) is valid if either 

p < q + 1, or 
p = q + l, ?7 > 1 



,bq;x) 
(3.4) 



(3.5a) 
(3.5b) 



are satisfied. It follows that pWq{x) exists for all GHS 
on the plane, but only for those with ry > 1 on the disk. 
The solution found, however, need not be positive. The 
positivity of pWq{x) must be checked in each particular 
case and may result in additional constraints on the pa- 
rameters. As was the case for the states, if I numerator 
parameters coalesce with I denominator parameters, the 
weight function pWq (x) reduces to the lower-order weight 
function p^iWq-i{x) of the remaining parameters, due to 
a well known property of the pFq and G functions. On 
the other hand, this can be seen also as a consequence of 
the defining condition H3.2|l . 

The resolution of unity, Eq. (|3.1|) . can be used to in- 
troduce a new representation by sandwiching it between 
two arbitrary state vectors, ($| and l^*), and writing the 
resulting scalar product in the form 



($1*) = i 1 d^Zp<^l{z)p^q{z), 



(3.6) 



with the wave functions in the GHCS basis defined ac- 
cording to 



("l*>- (3.7) 



We call this the GHCS representation of the state \^). 
The scries appearing in Eq. (|3.7|) defines an entire an- 
alytic function of C = -z*, which may be regarded as 



another representation. We call it the generalized hy- 
pergeometric analytic representation of l^*): 



p*,(c) = E 



.=0 VpPqH 

The scalar product (|3.t)|) then writes as 

($1*) = J pd^lq{0 p^iO p^qiO, 



(3.3) with the measure given by 



odfIq{Q = 



(3.8) 



(3.9) 



(3.10) 



Note that the analytic representations and the measures 
depend explicitly also on the parameter sets (ai, . . . , ap) 
and (6i,...,6g). For a given state we then obtain 
infinitely many analytic representations on the plane and 
unit disk bccing, respectively, elements of Bargmann and 
Hardy spaces with measure (|3.1UII . The usual Bargmann 
representation based on the coherent states is recovered 
for p = (jr = 0. The simplest Hardy space representation 
is based on the coherent phase states given in Eq. (|3.20() 
below II]. 

Next we consider the GHS on the unit circle and write 
z = e^'^ . We distinguish two types of states, denoted by 
]p; q] e*"^) and ]p; q; Lp) depending on whether the conver- 
gence criterion H2.5cll is satisfied or not. The normalized 
states [r] < 0) are defined by 



|p;<z;c'^^) 



^ oc 

—y 



\n) (3.11) 



and the unnormalizable states (rj > 0) by 



^ oo 



\n). 



(3.12) 



The normalization function in Eq. H3.11|l is a constant 
given by the generalized hypcrgcomctric function of unit 
argument. These states do not yield a resolution of unity, 
as the resulting Hausdorff moment problem 



27r 



dip pWq(ip) 



,2(71 — 771 



pPq{n) Sn 



(3.13) 



for the sought for function pWqlip) has no solution: 
Eq. H3.13|l implies the vanishing of all Fourier compo- 
nents of the function pWq{Lp) for n ^ m, thus reducing it 
to a constant which cannot match the values on the r.h.s. 
for all n = m. The lowest possible states ]l;0;c"^) ~ 
]a; • ; e*'^)(i.o) are not allowed, as the only parameter a 
cannot satisfy the conditions (|2.5c() and H2.6|l simulta- 
neously. The lowest normalizable states are ]2; l;e*''') = 
]ai, 02; b] e^^){2a) for ^1 = ai + a2 — b < in addition to the 
constraints (|2.6|) ; their normalization constant is given by 



4 



the Gaussian hypergeometric function of unit argument, 
(ai, a2; b; 1) = r(6)r(6- ai - a2)/[r(6- ai)r(&- 02)]- 
On the other hand, the simplest unnormaUzable states 



are 



|1;0;^) 



i;0) 



which in the hmit 



the well known phase states 



1 yield 



(3.14) 



T1=0 



These states, corresponding to ipoia; ■ ; n) |a— >i= 1 for all 
n, are the only ones yielding a resolution of unity with 
1^0 = l/(27r) according to Eq. (|3.13|) . 

We conclude this section by giving examples of GHCS 
for particular values of p and q. We quote for each state 
the parameter, the normalization and the weight func- 
tions. The complex variable z is denoted by a, \a\ < 00, 
for the GHCS on the plane {p < q + 1) and by e, |e| < 1, 
for the GHCS on the unit disk {p = q + 1). 

(a) The conventional coherent states |a) = |0;0;q!) 
correspond to p = q = 0, in which case oPoin) = 
oi^o(C) = and oWo{\a\^) = 1. 

(b) The states |0;l;a) = |- ; a)(0:i)i related with 
oPi{-;b;n) ~ n\(b)n for 6 > due to Eq. (|2.6|) . have 
the normalization function 

oFi(.;&;C) = r(6) /6-i(2yC) (3.15) 

and the weig ht function (0, p. 196, formula (5.39)]) 

0Wi{-;b;\af ) = 2 h-ii2\a\) Xb_i(2|a|), (3.16) 

where I^{x) and K^{x) are the modified Bessel functions. 
Since they are real and positive when v > —\ and x > 
[To| . the weight function o^i is positive for 6 > 0. 

(c) For the states |l;l;a) = |a;6;a)(i;i) with 
ipi(a; b; n) = n!(6)„/(a)„, the positivity constraints (|2.6(l 
require a and b to be either both positive or both nega- 
tive with equal negative integer parts. The normalization 
function is given by Kummer's confluent hypergeometric 
function iFi(a;b;C,) and the weight function by ([^ p. 
716, formula (9)]) 



iwi{a;b] \a\^) 



r(a) 

m 



,Fi{a-b-M^) X 



e-l"l' *(a-6;2-&;|an, (3.17) 



with 5'(a;&;C) the Tricomi confluent hypergeometric 
function. For a = b the above results reduce to those 
for the coherent states |0;0;q!) in example (a). 

(d) The states |l;0;e) = |a; • ; e)(]^.Q) originate from 
iPo{a; ■;"-)= nl/{a)n, with a > due to Eq. H2.6|l . The 
normalization function is given by 

iFo(a;-;C) = (1 - 0""^ (3.18) 
and the weight function by (0, p. 195, formula (5.35)]) 

a- 1 



iWo(a; • ; |e| ) 



ii-\e\^y 



(3.19) 



with a > 1 due to Eq. Ij3.5b|) . in which case also iWq > 0. 

The states |a; • ; e)(i;i) for a = 2k and Bargmann index 
k = i,l,|,... yield the Perelomov coherent states of 
the discrete series of the SU(1,1) representations, and for 
arbitrary non-negative k those of its universal covering 
group 0, 0] . Of particular interest are the states for 
fc = i, also known as coherent phase states ^3 



a: • ;e 



(i;0) 



(3.20) 



n=0 



Their weight function H3.19|l vanishes and the resolution 
of unity is applicable only as a limiting equation |l4j . 

(e) For the states |2;l;e) = |ai, 02; 6; e)(2;i) with 
2Pi( ai, 02; b; n) = n!(&)n/[(ai)„(a2)„], the constraints 
(|2.6() allow not only positive parameter values, but also 
pairs of negative-valued parameters with equal negative 
integer parts as well as a complex- valued pair (02 = a\). 
The normalization function is given by the Gaussian 
hypergeometric function 2^1(0-11 '3'2; C) ^-nd the weight 
fimction by (0, p. 198, formula (5.50)]) 



2^1(01,02; 6; |e| 



r(ai) r(a2) 



r(6) r(ai + a2-6-l) 



(l~kn°^+"=-''-'2Fi(ai,a2;6;|enx 
2Fi(a2 - 6, ai - 6; ai -f a2 - 6 - 1; 1 - jep), (3.21) 



with fli -|- 02 — 6 > 1 due to Eq. (|3.5b|l . It is not clear if 
this condition together with the constraints (|2.6|l are suf- 
ficient to guarantee the positivity of 2^1 in all cases. For 
coalescing parameters (oi = b or 02 — b) we recover the 
results of the previous example (d) for the states |1; 0; e). 



4. GENERALIZED HYPERGEOMETRIC 
STATES AS EIGENSTATES OF LOWERING 
OPERATORS 



In this section we introduce ladder operators of raising 
and lowering type, the latter having as eigenstates the 
generalized hypergeometric states introduced in Section 
121 These operators generalize the creation and annihila- 
tion operators of the harmonic oscillator in the case of 
the GH(C)S on the plane, and the Susskind-Glogower 
exponential phase operators in the case of the GH(C)S 
on the unit disk and circle. 

We define generalized hypergeometric lowering and 
raising operators by 



A = E p/^W \n){n+ll (4.1a) 

00 

vUl = E l"+l)H, (4.1b) 



n=0 
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where the coefficients pfq{n) are given by 

pfqi'iT') = pfq{ai,...,ap;bi,...,hq;n) 



.L+i) ;"+^\)---i ^±M, (4.2) 



with afa{Ti) = \/n + 1 and all products in the numerator 
(denominator) replaced by 1 if g = (p = 0). It is very 
satisfactory that these coefficients are real and positive 
due to the same constraints (|2.6|l already observed. The 
ladder operators depend explicitly on the two parame- 
ter sets, e.g. pUq — pUq{ai, . . . ,ap]bi, . . . ,bq). Due to 
the product structure of pfq{n), if I numerator param- 
eters coalesce with I denominator parameters, the oper- 
ator pUq reduces to the lower-order operator p-iUq^i of 
the remaining parameters. The ladder operators obey 
the noncanonical commutation relation 

oo 

[pUq, pU^] = (f/?H'-p/9("-l)') (4-3) 

n=0 

where p/q(— 1) = by definition, but also by Eq. H4.2|l . 

The action of the ladder operators on the Fock basis is 
given by 

pUq\n) = pfq{n-l) \n-l), (4.4a) 
pUl \n) = pfqin) (4.4b) 

justifying their name as lowering and raising operators, 
respectively. It is then easy to show that the GH(C)S are 
eigenstates of the lowering operator 



pUq \p;q;z) = z \p;q]z), 



(4.5) 



with complex eigenvalue z. Here the parameter functions 
pPq{n) defining the eigenstates and the coefficients pfq{n) 
defining the lowering operator are related by 

pPqin) = [pfqi0)pfq{l)pfqi2)--- pfqin-l))'. (4.6) 

Seen from this perspective the generalized hypcrgcomct- 
ric (coherent) states are also instancies of the very large 
class of so-called nonlinear coherent (or f-coherent) states 
0, 0, 01 . These are eigenstates of lowering operators 
of the type (|4.1a|) with pfq{n) replaced by ^/n + 1 /(n), 
where the f-parameters /(n) are quite arbitrary. In our 
case f{n) is given by Eq. H4.2|l without the factor v^rT+T. 



For p < 5 -I- 1 we denote pAq 



pUq and pA\ 



pUq, which may be considered as a generalization of the 
usual annihilation (a) and creation (a^) operators of the 
harmonic oscillator. In particular, we have 



0^0 



n+1 \n){n + 1| 



(4.7) 



It is then natural to call pAq (pAj) the generalized hyper- 
geometric annihilation (creation) operators. The eigen- 
states oi pAq are the GII(C)S on the plane, \p; q; a), with 



= pUq and pE}^ 



complex eigenvalue a, \a\ < oo. They may be considered 
as generalizations of the conventional coherent states |a), 
which are eigenstates of a. 
For p = q + 1 we write p 

The eigenstates of pEq are: (a) the GH(C)S on the unit 
disk, \p;q;e), with complex eigenvalue e, |e| < 1; (b) the 
GHS on the unit circle with eigenvalue e'"^ for both, the 
normalized states \p;q;e^'^), Eq. and the unnor- 

malizable states \p\ g; ip). Eg. (|3.12|) . This is similar to 
the Susskind-Glogower |l7l | exponential phase operator 
E, having as eigenstates the coherent phase states (I3.20f) 
on the unit disk and the phase states H3.14|l on the unit 
circle. Now the operator E is a limiting case of the op- 
erator iEqIu) for a — *■ 1: 



iEo{a) \a- 



E 



n=0 



\n) {n + l\ 



(4.8) 



This correlates with the limit for the corresponding eigen- 
states in Eqs. (|3.20|l and H3.14|l . It is then justified to 
consider the operators pEq and their eigenstates, the 
GHS on the unit disk and unit circle, as generalizations of 
the Susskind-Glogower exponential phase operator E and 
its eigenstates, the coherent phase states and the phase 
states, respectively. For this reason we call pEq and pE'^ 
the generalized hypergeometric exponential phase oper- 
ators. 

Given the nonhermitian ladder operators we may de- 
fine hermitian combinations by pQq = (pAj -I- pAq)/\/2 

and pPq — i{pA'^g — pAq)/^/2, and similarly by pCq = 

{pEl + pEq)/2 and pSq = i{pEl - pEq)/2. Our no- 

tation suggests to "interpret" pQq {pPq) as generalized 
hypergeometric coordinate (momentum) operators, and 
pCq ipSq) as generalized hypergeometric cosine (sine) op- 
erators (a la Susskind-Glogower). However, these opera- 
tors do not satisfy the usual commutation relations and 
we have not verified that they have the right spectrum. 
Their study is beyond the scope of this paper. 



5. PHOTON NUMBER STATISTICS 

The generalized hypergeometric (coherent) states lead 
to classical as well as nonclassical behaviour in the photon 
number statistics. The considerations in this section do 
not depend on the existence of a weight function and 
apply, therefore, to both GHS and GHCS. 

The GH(C)S defined in Eq. (EU have the Fock repre- 
sentation 



{n\p;q]z) = 



y/pPq{n) pA/'gdzP) ' 



(5.1) 



from which the photon number distribution follows 
1 



P\p;q-z)(.n) 



pPq{n) pNq{x) 



{x=\z\'). (5.2) 
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The factorial moments n^'^^ = {n{n — 1) • • • (n — fc + 1)) 
are then easily calculated in terms of the derivatives of 
the normalization function 



0.4- 



P"'q 



ix) 



pNq {x) da;*-' 



(5.3) 



In our case pMq = pFq and the derivatives of pFq can be 
expressed in terms of the same pFq with shifted parame- 
ter values un . We then obtain for the factorial moments 



pU^ '^ (x) — 
pFq{ai + k, 



X ~~7Z T 7~ T X 



(Qi)a 

• ■ ■ ibq)k 

^Qp + k;bi + k, . . . ,bq + k;x) 



pFq{ai, . . . ,ap;bi, . . .,bq;x) 



The mean photon mmrber, n = n^'^\^ is given by 

_ _ ai-'-Qp p-F,(ai + l,...,ap + l;bi-|-l,...,bg + l;a:) 



(5.4) 



bi---b„ 



oFg(ai,...,aj,-bi,...,bg-x) 



(5.5) 



and the Mandel parameter, Q ; 



,(2) 



/n[Tl[Ti,by 



ai---ap pFg{ai + l,...,ap + l-,bi + l,...,bg+l;x) 



fci ■■■()„ 



oFg{ai,...,ap:bi 



,bg-x) 



(6i + l)---(b, + l) pFg(ai + l,...,ap + l;bi+l,...,bg + l;x) 



(5.6) 



Positive, vanishing, or negative values of the Mandel pa- 
rameter correspond to super-Poissonian, Poissonian, or 
sub-Poissonian statistics, respectively. The latter is a sig- 
nature of nonclassical behaviour. If in the above expres- 
sions we let I a-parameters coalesce with I ^-parameters, 
we obtain the corresponding expressions for the lower- 
order states \p — l',q~ l\z). 

We now consider the photon number statistics of some 
particular GHCS. For each state we quote the photon 
number distribution, the mean photon number, and the 
Mandel parameter. We also show and discuss graphical 
representations for selected values of the parameters and 
compare them with those for the coherent states. 

(a) We start with the coherent states \a). Their photon 
number distribution is the Poisson distribution P\a) (n) ~ 

(|ap"/ri!) e~l"l', yielding oUq — \a\'^ and oQo — 0. 

(b) For the GHCS |0;l;a) : 



; b; a)(o;i) we have: 



^0;l;a)("-) 



a 



2n+b-l 



0"! 



n! r{b + n) /b-i(2|Q|)' 



Ib{2\a\) 



'6+1 



Ib-ii2\a\y 



(2|«|) Ib{2\a\) 



Ibi2\a\) Ib-ii2\a\) 



(5.7) 



(5.8) 



(5.9) 



The photon number distributions of the states |0;1;q;) 
in Fig. ^ peak at low values of the photon number n, 
their mean photon numbers in Fig. |21 increase for large 
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FIG. 1: Photon number distributions of the GHCS |0; 1; a) 

\- -jb; ce)(o;i) and the CS \a) for |a| = 3. 




a 



FIG. 2: Mean photon numbers of the GHCS |0;l;a> 
|- ;6;Qf)(0;i) and the CS \a) as a function of \a\. 




\a\ 



FIG. 3: Mandel parameters of the GHCS |0;l;a> 
|- ;6;a)(0;i) and the CS \a) as a function of |q|. 



7 




enough |a| much slower than for the coherent state, and 
the negative values of their Mandel parameters in Fig. O 
indicate nonclassical behaviour for all lal. 



(c) For the GHCS |1;1; 



we find 



-P|l;l;a>('^) 



(«)n 

n\ (&)„ iFi(a;&;|a|2)' 



ini 



a ■ 



1?7.1 



iFi(a + l;5+l; 



iFi(a;&;|a|2) 



j+l ifi(a+2;fc+2;|a|^) 
b+l iFi(a+l;6+l;|Q|^) ' 



(5.10) 



(5.11) 



(5.12) 



These relations reduce for a = 6 to those of the coherent 
states. This remark is very helpful in understanding the 
behaviour of the photon number distributions in Fig. 0] 
the mean photon numbers in Fig. |S1 and the Mandel pa- 
rameters in Fig. El All curves coincide for a ~ b with the 
CS-curves and depart from them as soon as a becomes 
different from b. The deviations are on cither side with 
respect to the CS-curves, depending on wether a < 6 or 
a > b. As an example, the Mandel parameters in Fig. El 
are positive (negative) for a < 6 (a > 6). 



(d) For the GHCS |l;0;e) 



-P|1;0;£>('^) 



(a), 



(1 



/(i;0) 



we have 



(5.13) 




-0.4 ■ 



3 

lal 



FIG. 6: Mandel parameters of the GHCS |l;l;a> 
|a;6;Q)(i.i) and the GS |a) as a function of |q|. 



states for this reason [lj|. The photon number dis- 
tributions arc shown in Fig. [3 and compared with the 
Poisson distribution for the same value of the variable 
(|e| = I = |a|). The mean photon number and the Man- 
del parameter have the same simple dependence on |e|. 
The Mandel parameter is independent of the parameter 
a and always positive (classical behaviour). 

(e) For the GHCS |2;l;e) = |ai, 02; e)(2;i) we obtain 



a iQo = 



1 



(5.14) 



These expressions make sens for any a > 0. Putting 
a = 2k and k = i,l,|,... in Eq. H5.13|l . we obtain 
the negative binomial distribution of the corresponding 
SU{1, 1) coherent states, also called negative binomial 



P\2;l;e)in) 



(ai)n (a2)r 



|2n 



7i! {b)n 2Fi{ai,a2;b; |eP) ' 



(5.15) 



. .2 aia2 2Fi{ai + l,a2 + l;b+l;\e\^) 
b 2i^i(ai,a2;0; e^) 
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n 

FIG. 7: Photon number distributions of the GHCS |1; 0; e) = 
|a; • ; e)(i;o) and the CS \a) for |e| = | = 



0.5 




-3-2-10 1 2 3 

e-cp 

FIG. 8: Husimi phase distributions P^^l^^ (6) of the GHCS 
|0;l;a) = h&;a>(0;i) and of the CS |q) for \a\ = |. 



2Q1 



|2 (ai + l)(a2 + l) 2fl(ai+2,a2+2;b+2;|er) 
I (6+1) 



2Fi(ai + l,a2+l;6+l;|e|2) 

(5.17) 

For b = a2 (or 6 = ai) we recover the results of the 
previous example (d). We refrain from grapliical repre- 
sentations as there are too many parameters. 



6. HUSIMI AND HUSIMI PHASE 
DISTRIBUTIONS 

In this section we consider the Husimi distributions 
(Q functions) and the corresponding phase distributions 
for the generahzed hypergeometric (coherent) states. We 
also comment on the Pegg-Barnett phase distributions. 

The Husimi distribution for the GH(C)S \p; q; z) is 
given by the modulus squared of its wave function in 
the coherent state basis la) 



Q\p;q:z){a) = - | (a |p; (?; z) | ^ , 



where 

{a\p;q;z) 



E 



(6.1) 



(6.2) 



It is positive, normalized to 1 with measure (Pa and pro- 
vides a two-dimensional probability distribution over the 
complex a-plane. Writing a — |a|e*^ and integrating 
over the radial part \a\, the Husimi phase distribution is 
obtained in the form I2II 



n.n'—O 

g^^\n,n') e-<^-'^')\ (6.3) 



where the Fock components are given by Eq. (|5.1|l and 
the newly introduced coefficients by 



r( 



n-[-n 



1) 



^T{n + l) r(n' + l) 



(6.4) 



In a different approach the phase states (|3.14|l can be 
used to define the phase distribution 



{PB) 



{e) = mp-q-z) 



(6.5) 



It has the same form as Eq. (|6.3(l . with Q^'^\n^n') re- 
placed by g^^^\n,n') = 1. Since the same result is ob- 
tained for the (differently motivated) Pegg-Barnett phase 
distribution |22l |. we have used the upper index (PB) in 
Eq. H6.5|) . Note that the dependence of the phase dis- 
tributions on the phases is through the phase difference 
(9 — if) only, where if is the phase of the complex variable 
{z = Izle**^) defining the GH(C)S. 

We display in Figures ISlllOl the Husimi phase distribu- 
tions of the GHCS \p;q;z) for {p;q) — (0;1),(1;1), and 
(1;0), respectively, and for different values of the corre- 
sponding parameters. All distributions peak at 9 = 
The heights of the peaks increase: with decreasing b in 
Fig. IHl with decreasing b and increasing a in Fig. El and 
with increasing a in Fig. EH In Fig- IHI the phase dis- 



tributions Vl^l^^{e) for a = 5 and Vl^^'{9) would be 
identical; for a > 6 (a < b) the heights of the peaks then 
move above (below) the CS peak. This is similar to the 
behaviour observed in the photon number distributions 
of the states |1; 1; a). 

The Pegg-Barnett phase distributions look qualita- 
tively similar, so we do not display them explicitly. The 
peaks at 9 = ip are always higher (due to C/^^^^(n, n') > 
Q^'^\n^n')) and consequently the distributions are nar- 
rower in order to preserve the normalization. 
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Husimi phase distributions Vy^\.^^{B) of tlie GHCS 
= |a;6;a)(i.i) and V'^^^iQ) of the CS \a) for |q| = |. 
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FIG. 10: Husimi phase distributions V\2l.^-. (0) of the GHCS 



|l;0;e> = |a;.;e)(i;0) and pff^'(e) of the CS |a) for |e| = f = 



7. GENERALIZED HYPERGEOMETRIC 
HUSIMI AND GENERALIZED 
HYPERGEOMETRIC HUSIMI PHASE 
DISTRIBUTIONS 



By analogy with the conventional Husimi distribution 
we define in this section the generalized hypergeomet- 
ric Husimi distribution of a normalized state by the 
modulus squared of its wave function in the GHCS basis, 
Eq. JTtII : 



Q[;f(z) ^ Q[^f (ai,...,ap;6i,. 



.,bq;z) 



pWgdzp) {p;q;z\'4>){'ip\p;q;z). (7.1) 



It provides a two-dimensional true probability distribu- 
tion over the complex z-plane for p < q + 1 and over the 
unit z-disk for p = q + I. The corresponding normaliza- 
tion condition 



= 1 



(7.2) 



is a consequence of the resolution of unity, Eq. H3.1|l . It is 
therefore mandatory for the state \p; q; z) in the definition 
(|7Hl to be a GHCS; the state \'\\)) can be any normalized 
state, in particular a GHS or GHCS. For later reference, 
we call the states j?/') and \p\ q; z) the signal and analyzing 
states, respectively. In defining generalized hypergeomet- 
ric Husimi distributions we have formally replaced the 
analyzing CS in the usual Husimi distribution by an an- 
alyzing GHCS and have included the weight function into 
the distribution rather than into the integration measure. 
The conventional Husimi distribution is then recovered 
for p = (? = 0, i.e. Q\^){a) 



The generalized hypergeometric Husimi distributions 
can be used to define corresponding phase distributions 
by integrating over the modulus of the complex variable 

z = |z|e*« {x = |z|2): 



(7.3) 



where i? = oo (i? = 1) for the GHCS on the plane (disk). 
Performing the integral with the help of Eq. (|3.2|) , obvi- 
ously extendable to arbitrary (here halfinteger) n, the 
generalized hypergeometric Husimi phase distributions 
are obtained in the form 



-. oc 



i{n-~n )6 



They are normalized according to 



(7.4) 



(7.5) 



where the integration is over the 27r-interval chosen for 
the phase 6*, usually [— tt, tt] or [0, 27r]. 

The phase distributions H7.4|l and (|6.3() have the 
same formal structure, only the coefficients are different. 
Those appearing in Eq. (|7.4|) are given by 



(71, n') 



pPq{ 2 ) 



\/pPq{'n) pPq{n') 



(7.6) 



in terms of the moments pPq{n) from Eq. (|2.2() . They are 
symmetric in n and n' and normalized to 1 for n = n' , 
thus insuring the normalization condition l|7.5|l . These 
coefficients encode the information on the analyzing state 
\p]q;z). The information on the signal state \ip) is con- 
tained solely in the product ipn 4^^', where ^pn = {n-W- If 
the system is described by a density operator p, this prod- 
uct should be replaced by the density matrix elements 
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in the Fock basis, pn,n' = {nlpW)') correspondingly, p 
should appear instead of in Eq. H7.1|l . If the de- 

pendence on the parameters is needed, we shall write 



lOi, . . . , Op; oi, . . . , 

and ^^(P'«)(ai, . . . ,ap;6i. 
The latter can be written also in factorized form 



) for the phase distributions 
bq] n, n') for the coefficients. 



t(p;9) 



{n,n') 



with the one-parameter coefficients given by 

r(^ + a) 



y^r{n + a) r{n' + a) 



(7.7) 
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Comparing the product structures of pPq{n) and 
Eqs. and (jZ2I), it follows that 

each Pochhammer symbol (a)„ contributes a coefficient 
^^("'(71,71'); the factor Q'-'^^ {n,n') = Q^^'>{n,n') comes 
from T{n + 1) = (1)„ according to the same rule. If 
in Eq. 1)7.7(1 p or g or both are zero, the corresponding 
product(s) should be replaced by 1. Thus C/^"-") (n, n') = 
that the conventional Husimi phase dis- 
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FIG. 11; Generalized hypergeometric Husimi phase distribu- 



tions Vi 



the CS ja) for |a| 



0.6- 



; b; 9) and Husimi phase distribution V^^^ (9) of 



tribution 



IS given by Pl^^\e) = Vl^'{0). Simi- 



,(0;0), 



larly, the Pegg-Barnett phase distribution is given by 



(0) 



p(i;o)/ 



since C/^^-") (a; • 



n, n 



Q^'^\n,n')/Q^"-^ {n,n') equals 1 for a = 1. More generally, 
the generalized hypergeometric Husimi and correspond- 
ing phase distributions of order {p — I; q — I) result from 
those of order (p; q) by letting I numerator parameters 
coalesce with / denominator parameters. 

We now consider some generalized hypergeometric 
Husimi phase distributions for the Fock and coherent 
states. For any Fock state |A^) we have ipn = 6n.N and 
the phase distributions l|7.4|) become uniform with value 
l/(27r), as expected. 

For the coherent states |a) with a = |a|e"^ and 
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l"l V2 |Q,|n Q^nv ^^6 rcsuhing phasc distribu- FIG. 12: Generalized hypergeometric Husimi phase distribu- 



tions ((7.4(1 depend only on the phase difference {Q — Lp). 
We display in Figures I11II13I the generalized hypergeo- 
metric Husimi phase distributions of a coherent state for 
(p;^) = (0; 1), (1; 1), and (1;0), respectively. All distri- 
butions peak at the coherent phase 9 = Lp. The heights 
of the peaks decrease: with decreasing h in Fig. ^] with 
decreasing h and increasing a in Fig. 1121 and with in- 
creasing a in Fig. In Fig. the phase distribution 
^|q)^' (a; b] 9) for a ^ b and the Husimi phase distribution 
"^{a)^^) would yield the same results; for a < 6 (a > b) 
the heights of the peaks then move above (below) the 
peak of 'P^^^\9). This behaviour is similar to the one 
observed in the conventional Husimi phase distributions 
of the states |1; 1;q!) as well as in their photon number 
distributions. 

The phase distributions shown in Figures I11II13I are 
"dual" to those in Figures ISlH 01 in the sense that the sig- 
nal and analysing states have been interchanged. The 



tions Pi^^J^' (a; 6; and Husimi phase distribution V^^) (9) of 
the CS \a} for \a\ = |. 



signal (analyzing) state is a GHCS (CS) in Figures ISt llOl 
and a CS (GHCS) in Figures E] [131 In the first case we 
deal with the usual Husimi phase distribution V 



(0;0) 



of 

the GHCS \p;q;z), in the second case with the general- 
ized hypergeometric Husimi phase distribution 7^^^!^^"^ of 
the conventional CS \a). The two sets of figures show an 
opposite behaviour in the sense that for the same varia- 
tion of the parameters the heights of the peaks increase 
in one set and decrease in the other set. 

The phase distributions showing this dual behaviour 
originate from generalized hypergeometric Husimi dis- 
tributions of the type (5[^.'|.£) (2^), where both the 
signal {\p;q;z)) and the analysing {\p;q;z)) states 
are GHCS and therefore interchangeable. We call 
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FIG. 13: Generalized hypergeometric Husimi phase distribu- 
tions P'^j"' (a; • \ 0) and Husimi phase distribution P^^-^\0) of 
the CS |q) for \a\ = |. 



the generalized 



hypergeometric Husimi distribution 
bi, . . . ,bq] z) and the one with intcr- 

6„; z), dual to 



changed states, Qip.'^?^) (ai, . . . , a^; foi, . . . , 6,; z), 
each other and extend this definition to the correspond- 
ing phase distributions, v\^:^\.^{e) and v\^:l\^{e). Dual 
distributions arc nontrivial due to the nonorthogonality 
of the GHCS and this is true even if for p = p and q ~ q 
the two sets of tilded and untilded parameters are identi- 
cal. In this case the self-dual distributions Q\p.q\) (z) are 
determined by the scalar product in Eq. 12. 4|) . i.e. 



1 pWq{\z\^) \yMq{z*~z)f 

TT ,M,{\Z\^) ,M,i\~zY) 



(7.9) 



Self-dual distributions actually provide a measure of the 
nonorthogonality of the GHCS with identical parame- 
ter sets. The simplest self-dual distribution is the usual 
Husimi distribution for a conventional coherent state, 
given by the well known overlap of the two coherent 
states: = :^|(Q:|a)P = Q\a)[pL). The correspond- 

ing self-dual phase distribution is (&) , also shown in 
Figures ISIIT^ for comparison. 



8. SUMMARY 



In this paper we have introduced a large class of 



holomorphic quantum states |ai. 



,ap; 6,; z), de- 



pending on a complex variable z and two sets of parame- 
ters, (oi, . . . , Op) and . . . , 6,). We called them gener- 
alized hypergeometric states (GHS) as their normaliza- 
tion function is given by the generalized hypergeometric 



; 6i, . . . , 6g; |z| ). Depending on 



function p_F'q(ai,.., 
the domain of convergence of the latter, the GHS are 
defined on the whole z-planc (|z| < oo), the open unit 
disk {\z\ < 1), or the unit circle (|z| = 1). Accord- 
ingly, these states may be considered as generalizations 
of the usual coherent states, the coherent phase states, 
and the phase states. Those GHS yielding a resolution of 
unity with a positive weight function define the subclass 
of generalized hypergeometric coherent states (GHCS). 
We have determined the corresponding weight functions 
by Mellin transform techniques, but their positivity must 
be checked separately. The GHS on the plane (unit disk 
and circle) are eigenstates with eigenvalue z of suitably 
defined lowering operators pAq (pEq); the latter may be 
considered as generalizations of a (E), the usual annihi- 
lation operator of the harmonic oscillator (the Susskind- 
Glogower exponential phase operator). 

To reveal the physical content of the GH(C)S wc have 
studied their photon number statistics (photon number 
distribution, mean photon number, Mandel parameter) 
and phase properties (Husimi and Pegg-Barnett phase 
distributions). Numerical results were displayed for par- 
ticular states and parameter values. 

The GHCS on the plane (unit disk) can be used to de- 
fine new analytic representations of arbitrary quantum 
states, which are elements of corresponding Bargmann 
(Hardy) spaces. They can be used also to define general- 
ized hypergeometric Husimi functions over the complex 
plane (unit disk) and corresponding phase distributions. 
Of particular interest are the generalized hypergeomet- 
ric Husimi functions when considered for another GHCS, 
as then the two GHCS can be interchanged (dual distri- 
butions). The corresponding phase distributions show a 
correlated behaviour. Self-dual distributions are nontriv- 
ial and provide a measure of the nonorthogonality of the 
GHCS. 
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